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ABSTRACT 

The evolution of the galaxy stellar mass-star formation rate relationship (A/* — SFR) provides 
key constraints on the stellar mass assembly histories of galaxies. For star-forming galaxies, 
M*— SFR is observed to be fairly tight with a slope close to unity from z ^ ^ 2, and 
' (— I ' it evolves downwards roughly independently of M*. Simulations of galaxy formation repro- 

pi |. duce these trends, broadly independent of modeling details, owing to the generic dominance 

of smooth and steady cold accretion in these systems. In contrast, the observed amplitude of 
the 7\f*— SFR relation evolves markedly differently than in models, indicating either that stel- 
^ lar mass assembly is poorly understood or that observations have been misinterpreted. Stated 

' in terms of a star formation activity parameter asf = (M^/S'i^i?) /(inubbic — IGyr), models 

predict a constant agf ~ 1 out to redshifts z ^ 4+, while the observed il/*— SFR relation 
indicates that agf increases by ~ x3 from z ^ 2 until today. The low agf (i.e. rapid star for- 
mation) at high-z not only conflicts with models, but is also difficult to reconcile with other 
observations of high-z galaxies, such as the small scatter in A/*— SFR, the slow evolution of 
star forming galaxies at z ~ 2 — 4, and the modest passive fractions in mass-selected sam- 
I pies. Systematic biases could significantly affect measurements of AI^, and SFR, but detailed 

i considerations suggest that none are obvious candidates to reconcile the discrepancy. A spec- 

ulative solution is considered in which the stellar initial mass function (IMF) evolves towards 
more high-mass star formation at earlier epochs. Following Larson (1998), a model is inves- 
' tigated in which the characteristic mass M where the IMF turns over increases with redshift. 

I Population synthesis models are used to show that the observed and predicted Af*— SFR evo- 

lution may be brought into general agreement if A/ ~ 0.5(1 + z)^Mq out to z ~ 2. Such IMF 
evolution matches recent observations of cosmic stellar mas s growth, and the res ulting z = 
i cumulative IMF is similar to the "paunchy" IMF favored bv Fardal et al. I (12007') to reconcile 

' the observed cosmic star formation history with present-day fossil light measures. 
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1 INTRODUCTION 

How galaxies build up their stellar mass is a central question in 
galaxy formation. Within the broadly successful cold dark mat- 
ter scenario, gas is believed to accrete gravitationally into grow- 
ing dark matter halos, and then radiative processes enable the gas 
to decouple from non- baryonic matter and event ually settle into a 
star-forming disk (e.g. IMo, Mao, & Whit3ll998 '). However, many 
complications arise when testing this scenario against observations, 
as the stellar assembly of galaxies involves a host of other processes 
including star formation, kinetic and thermal feedback from vari- 
ous sources, and merger-induced activity, all of which are poorly 
understood in comparison to halo assembly. 

A key insight into gas accretion processes is the re cent recog- 
nition of the importance of "cold mode" accreti on jKatz et aU 
l2003l : lBimboim & Dekel1l2003l : lKeres et al. II2005I) . where gas in- 



falling from the intergalactic medium (IGM) does not pass through 
an accretion shock on its way to forming stars. Simulations and 
analytic m odels show that c old mode dominates global accretion 
at z > 2 jKeres et al. Il2005 | ). and dominates accretion in all ha- 
los w ith masses < lO^^Mr:^ jBirnboim & Dekel l2003l : lKeres et al.1 
l2005l) . 

The central features of cold mode accretion are that it is (1) 
rapid (2) smooth, and (3) steady. It is rapid because it is limited 
by the free-fall time and not the cooling time; it is smooth because 
most of the accretion occurs in small lumps and not through ma- 
jor mergers jMurali et al.ll200ll : iKeres et al. II2005I : IGuo & White] 
2007); and it is steady because it is governed by the gravitational 
potential of the slowly growing halo. A consequence of cold mode 
accretion is that the star formation rate is a fairly steady function 
of tim e (e.g. iFinlator et al.ll2006l ; [pinlator. Dave & Oppenheimer I 
l2007h . which results in a tight relationship between the stellar 
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mass Mt and the star formation rate SFR. Hence simulations of 
star-forming galaxies generically predict a tight relationship with 
M, ocSFR that evolves slowly with redshift. In detail, a slope 
slightly below unity occurs owing to the growth of hot halos around 
higher-mass galaxies that retards accretion. 

The stellar mass-star formation rate (Af*— SFR) relation for 
star forming galaxies has now been observed out to 2: '--^ 2, thanks 
to improving multiwavelength surveys. As expected from models, 
the relationship is seen to be fairly tight from z ~ 0—2, with a slope 
just below unity and a scatter of < .3 dex jNoeske et al. 112007a : 
lElbaz et al. IlioOTi : lOaddi et al. IboOTh . It also evolves downwards 
in amplitude to lower redshift in lock step fashion, suggesting a 
quiescent global quenching mechanism such as a lowering of the 
ambient cosmic density, again as expected in mod els. The range 
of data used to quantify this evolution is impressive: INoeske et al. I 
( l2007ij) used the AEGIS multi-wavelength survey in the Extended 
Groth Strip to quantify the A/* -SFR relation from z ~ 1 — » 0; 
lElbaz et al. I ilOOt u sed GOODS at z ~ 0.8 - 1.2 and SDSS 
spectra at 2 ~ 0; and lDaddi et al. I ( l2007h used star-forming BzK- 
selected galaxies with Spitzer, X-ray, and radio follow-up to study 
the relation sXz ~ 1.4 — 2.5. In each case, a careful accounting was 
done of both direct UV and re-radiated infrared photons to measure 
the total galaxy SFR, paying particular attention to contamination 
from active galactic nuclei (AGN) emission (discussed in [Q. Mul- 
tiwavelength data covering the rest-optical were employed to accu- 
rately estimate Af, . This broad agreement in A/, —SFR slope, scat- 
ter, and qualitative evolution between model predictions and these 
data lends support to the idea that cold mode accretion dominates 
in these galaxies. 

Yet all is not well for theory. Closer inspection reveals that 
the amplitude of the Af, —SFR relation evolves with time in a way 
that is inconsistent with model expectations. This disagreement is 
fairly generic, as shown in ij2l arising in both hydrodynamic simu- 
lations and semi-analytic models, and is fairly insensitive to feed- 
back implementation. The sense of the disagreement is that going 
to higher redshifts, the observed SFRs are higher, and/or stellar 
masses lower, than predicted in current models. 

The purpose of this paper is to understand the implications 
of the observed Af, — SFR relation for our theoretical view of how 
galaxies accumulate stellar mass. In ij3]it is argued that Af,— SFR 
amplitude evolution implies a typical galaxy star formation history 
that is difficult to reconcile with not only model expectations, but 
also other observations of high- 2: galaxies. Possible systematic ef- 
fects that may bias the estimation of SFR and Af, are considered 
in ^ and it is argued that none of them are obvious candidates 
to explain the discrepancy. Finally a speculative avenue for recon- 
ciliation is considered, namely that the stellar initial mass function 
(IMF) has a characteristic mass that evolves with redshift (ijS}. This 
model is constrained based on the Af,-SFR relation in fjS] The im- 
plications for such an evolving IMF are discussed in SjT] Results 
are summarized in ^ A ACDM cosmology with Q, — 0.25 and 
ffo ~ 70 km/s/Mpc is assumed throughout for computing cosmic 
timescales. 



2 Af, -SFR: MODELS VS. OBSERVATIONS 

To begin, the observed Af,— SFR relation is compared with 
model predictions, in order to understand the origin of the 
relation. The primary simulation s used are cos mological hy- 
drodynamic runs with Gadget-2 l lSDringel|[200^ incorporating 
momentum-driven outflows, which uniquely match a wide variety 



of IGM and galaxy observations ( 'e.g. Oppenheimer & Dave 1200^ : 



iDave. Finlator & Oppenheirner] l200fl : 
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2007h . 



The new simulations include stellar recycling assuming a Chabrier 
IMF, by returning gas (and metals) into the ambient reservoir, in 
addition to th e usual heating, cooling, and star formation, and feed- 
back recipes iSpringel & Hemquist II2003I : lOppenheimer & Dave I 
1 20061) . Outflow parameters are computed based on properties of 
galaxies identified with an on-the-fly group finder. These and other 
improvements that go towards ma king simulations more r ealisti c 
down to z = are described in lOppenheimer & Dave I ( l2007h . 
Two runs to z = are considered here, each with 256^ dark mat- 
ter and an equal number of gas particle, in cubic volumes of 32 
and 64/i~^Mpc (comoving) on a side. The assumed cosmology 
is ACDM with Qm = 0.3, = 0.7, ffo = 69 km/s/Mpc, 
fib = 0.048, and as — 0.83. The resolved galaxy population, 
defined here as those with stellar masses greater than 64 gas par- 
ticle masses, is limited by stellar masses of 2.5 X 10^ Af© and 
2 X lO^"Af0 in the two volumes (this is somewhat more conser- 
vative than in lFinlator et al. ||2006 . in order to ensure accurate star 
formation histories). 

To explore the impact of feedback parameters, simulations 
with two other outflow mode ls are considered, namely on e with 
no outflows, and one with the lSpringel & Hemquist I j2003h "con- 
stant wind" out flow model, taken from the suite of outflow 
runs described in lOppenheimer &Davel ( l200^ . These runs have 
32^~^Mpc box sizes with 2 x 256^ particles, and they were 
evolved only down to z = 2. Their cosmology is slightly differ- 
ent, having ag = 0.9 from first-year WMAP parameters, but this is 
not a large effect at the epochs considered here. 

In hydro simulations, star formation rates and stellar masses 
are obtained by summing over particles in galaxies identified us- 
ing the group finder SKID (Spline Kernel Interpolative DEN- 
MAX; http://www-hpcc.astro.washington. edu/tools/skid.html I; see 
iKeres et al. I l l2005 l) for details. Gas particles that have high enough 
density to form stars have instantaneous star formation rates com- 
puted in Gadget-2; it is this SFR that is summed to obtain the in- 
stantaneous SFR of a galaxy. Of course, observed galaxy SFR's and 
Af , 's are measured quite differently than in simulations, a point that 
will figure prominently in f|4] Note that these hydro simulations do 
not include any feedback from active galactic nuclei (AGN), or any 
process that explicitly truncates star formation in massive galax- 
ies, hence they do not reproduce the evolution of passive galaxies. 
Here the focus is on star-forming galaxies only, with the assump- 
tion that galaxies are quickly squelched to become passive and so 
do no t appear in star-forming selected samples (e.g. Salim et ah] 
l2007h . Since all galaxies in our simulations would be observation- 
ally classified 'as star-forming, no cuts are made to the simulated 
galaxy population. 

To complement the hydro simulations, results for Af, — SFR 
are taken from the Millenium simul ation plus semi-analyti c mode l 



Kitzbichler&Whitel ( l2007h . 



Elbaz et al. I ( l2007h . Besides 



("Milleni um SAM") Ughtcon es of 
as fit by iDaddi et al. I bOOTh and 
being a completely different modeling approach, the Millenium 
SAM also includes AGN feedback that broadly reproduces the ob- 
served passive galaxy population and its evolution. Additionally, 
a differ ent semi-analytic model ba sed on the Millenium simula- 
tion by iDe Lucia & Blaizoti i2007h was examined, to see if vari- 
ations in semi-analytic prescriptions can res ult in significant dif- 
ferenc es in Af, — SFR. It was found that the iDe Lucia & Blaizot I 
I2OO 7I) model produced a somewhat shallower slope than the 
Kitzb ichler & White I ^20071) model, which e.g. at 2 ~ 2 resulted 
in a smaller offset from observations at Af, ~ 10^" Af©, but a 
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Figure 1. M*— SFR relation for simulations versus observations at z = 
0,1,2. Red points show a 32/i~^Mpc hydro simulation box, magenta 
show 64/i~^Mpc box. A running mean (with standard deviations shown 
as vertical bars) is shown separately for each simulation; agreement in 
overlap region is excellent. Green lines show results fro m the Millenium 
SAM. Blue lines show fits to observations at 2: = 2 I D addi et al. 
and z = 0, 1 telbaz et al. 1 12007). with dashed lines indicating ±0.3 dex 
scatter jNoeske etam 2007a). The models' slopes and scatter are in good 
agreement, but the amplitude deviates from observations at high-2. Lower 
right: Star formation histories in 100 Myr bins for four star forming 
galaxies selected from the 32/i~^Mpc run closest to a final logM» = 
11, 10.5, 10, 9.5 (top to bottom curves). Note the lack of large bursts and 
the similarity of the star formation history shapes, which yields the predic- 
tion of a tight M, ocSFR (roughly) relation. The magenta line shows the 
form of the SFH assumed in PEGASE modehng discussed in ^ 



larger one at M, ~ 10 ^Mr. Since the re sults are qualitatively 
similar, only the lKitzbichler & White I i2007h models are discussed 
further. While there are many more semi-analytic models in the lit- 
erature, and perhaps even some that are in better agreement with the 
M*— SFR relation, it is beyond the scope of this work to consider 
them all. 

Figure [T] (first three panels) shows the M*— SFR relation 
at z = 0,1,2 in simulations, co mpared to pa rameterizations 
from obs ervations (bl u e line s) by lElbaz et alTI ([2007,) {z ~ 
0, 1) and iDaddietaJ.l (l2007h (z » 2). The observed scatter 
of 0.3 dex ('Noeske et al.ll2007al : lElbaz etal.ll200l lOaddi et al. I 
l2007h is indicated by the dashed blue lines. It is w orth noting that 
other z ~ 2 observations su ggest a larger scatter IShaplev et al. I 
I2OO5I : IPapovich et all I2OO6I) , so the observational picture on 
this is not entirely settled . Finally, the green lines show the 
iKitzbichler & White I ilOOH) Millenium SAM results, which gen- 
erally track the hydro simulation results, particularly in amplitude 
which will be the key aspect. 

At 2: = 0, the relations are generally in agreement, though 
the observed slope of M, ocSFR"'^^ is somewhat shallower than 
that predicted by the hydro simulations (~ 0.9). The Millenium 



SAM produces excellent agreement at z = 0, having broadly been 
tuned to do so. By z = 1, t he predicted hydro simulation slopes are 
similar to that observe d by Elbaz et al. '/2007) though steeper than 
seen bv lNoeske et al. 1 (2007a). but more importantly the amplitude 
is low by ~ x2 — 3. At 2 — 2, this trend continues, with a ~ 
x4 — 5 amplitude offset. The Millenium SAM has a slope that 
remains near the z = value, hence it is somewhat shallower than 
observed at 2: ~ 2, but the amplitude evolves slowly much like in 
the hydro simulations. The rapid evolution i n observed Af, — SFR 
amplitude has been n oted by several authors dPapovich et al. I2OO6I : 
iNoeske et aLll2007al) . a long with the f act that it is more rapid than 
in the Millenium SAM dElbaz et al. (2 007: Daddi et al. 2007). The 
new aspect demonstrated here is that the disagreement persists for 
current hydrodynamic simulations, and as will be shown in the next 
section this is true regardless of the main free parameter in such 
models, namely feedback implementation. 

The slope and scatter are a direct reflection of the similarity in 
the shape of galaxy star formation histories (SFH) among various 
models and at various masses. For a galaxy observed at time to, 
neglecting mass loss from stellar evolution. 



A/*(to)=/ SFR{t)dt = SFR{to) 



SFR(t) 
SFR(to) 



dt. (1) 



The integral depends purely on the form of the SFH up to tg. So 
long as the SFH form does not vary significantly with Al,, the 
slope of il/, — SFR will be near unity, and so long as the current 
SFR is similar to the re cent past-averaged one , the scatter will be 
small. As pointed out bv lNoeske etaTI ( i2007bh . the observed tight 
scatter of < 0.3 dex implies that the bulk of star-forming galaxies 
cannot be undergoing large bursts; this will be discussed further in 
Sj3] This is consistent with a fairly steady mode of galaxy assem- 
bly as generically predicted by cold mode accretion. Note however 
that the reverse conclusion does not necessarily hold, in the sense 
that a slope near unity is not uniquely indicative of a steady mode of 
galaxy forma tion; other scenarios ( such as the staged galaxy forma- 
tion model o f iNoeske et al. l2007bl) can be constructed to reproduce 
a slope near unity, as discussed in ^ 

Example galaxy SFHs from the 32/i~'^Mpc momentum- 
driven outflow hydro simulation are shown in the lower right panel 
of Figure [T] for a range of z = stellar masses from 10" M0 - 
lO^'^M©, in 100 Myr bins. They are remarkably self-similar, and 
generally quite smooth. The increased variance in smaller galaxies 
owes to discreteness effects in spawning star particles (these SFHs 
are computed from the history of galaxy star particles rather than 
the instantaneous SFR at each time); the true scatter in instanta- 
neous SFR increases only mildly at small masses, from ~ 0.2 dex 
at M* > 10^^ Mq to ~ 0.3 dex at M, > 10^'^ Mq, as shown in 
Figure [T] 

The early simulated SFHs are characterized by a rapid rise 
in star formation to z 6, corresponding to rapid early halo 
growth dLi et al. 1 1200% . Next, there is a period of roughly con- 
stant star formation from z ~ 6 — > 2; this will be the key 
epoch for the present work. Finally, at 2 < 2 the star forma- 
tion rates decline modestly, with more massive systems showing 
greater decline. This pattern reflects the changing a balance be- 
tween the growth of potential wells being able to attract more mat- 
ter, and cosmic expansion which lowers the ambient IGM den- 
sity and infall rates. It is qualitatively similar to the SFH inferred 
for the Milky Way. In the models, the fact that the SFH depends 
weakly on M* gives rise to M« ocSFR, and the fact that merger- 
driven st arbursts are not a significant growth pat h for galaxy stel- 
lar mass dKeres et al.1l2005l ; fGuo & White Il2007h gives rise to the 
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small scatter. In detail a slope slightly below unity occurs because 

more massive galaxies tend to have their accretion curtailed by 

f II 

hot halos (e.g. Keres et al. 2005), resulting in "natural downsiz- 
ing" [Neistein, va n den B osch & Dekel 2006). 

Note that the SFRs of individual galaxies in simulations do 
not drop by a large amount over a Hubble time. Even for the most 
massive star forming galaxies today, the SFR in simulations drops 
only be a factor of few since z = 2, and smaller galaxies show 
essentially no drop. Meanwhile, the amplitude of M,-SFR drops 
by nearly an order of magnitude in the models, and even more in the 
data. So the amplitude evolution owes predominantly to the growth 
of M, of individual galaxies. In other words, models indicate that 
the A/,— SFR relation doesn't evolve downwards in SFR so much 
as upwards in AI, ■ 

Overall, while the slope and scatter of the M, — SFR relation 
are in broad agreement, there is a stark and generic discrepancy in 
the amplitude evolution between the models and data. The agree- 
ment in slope and scatter is an encouraging sign that our under- 
standing of gas accretion processes is reasonably sound. Examining 
the discrepancy in amplitude evolution is the focus of the remainder 
of this paper. 



3 EVOLUTION OF THE STAR FORMATION ACTIVITY 
PARAMETER 

The amplitude evolution of A/* — SFR can be recast in terms of the 
timescale for star formation activity. Models generally predict that 
galaxies have A/*/SFR~ tn, where tn is the Hubble time, perhaps 
reduced by 1 Gyr because few stars were forming at 2 > 6, as seen 
e.g. in the lower right panel of Figure[T] 

A convenient parameterization is provided by the star forma- 
tion activity paramete^ defined here as the dimensionless quantity 
Qaf = (Af,/SFR) / (tH — IGyr). Physically, this may be regarded 
as the fraction of the Hubble time (minus a Gyr) that a galaxy needs 
to have formed stars at its current rate in order to produce its current 
stellar mass. Note that oisi is in general a function of Af, . 

Figure |2] shows Usiiz) for simulations and observations, top 
panel showing values at Af* — IO^^'^A/q, bottom showing A/, = 
10^" A/0. The values for models have been obtained by fitting 
power laws to simulated galaxies' A/*— SFR. The momentum- 
driven outflow hydro simulation results are shown as the red cir- 
cles, while cyan an d magenta circles show results for no winds and 
the wind model of Springel & Hemquist ( 2003) at z = 2 — 4, 
respectively. The outflow model has only a mild effect on Osf, 
because the form of the SFHs are broadly similar regardless of 
outflows, even though the total amount of stars formed is signifi- 
cantly different (Oppenheimer & Dave 2006). The most dramatic 
outlier among models is the high-mass bin of the no-wind simu- 
lation, where the rapid gas consumption in large galaxies results 
in lower current star formation rates relative to the past-averaged 
one (i.e. natural downsizing); this actually makes Qgf larger and the 
discrepancy with observations worse. Millenium SAM results are 
quite similar to the hydro runs, again showing little evolution in 
Osf , although the shallower slope results in a somewhat smaller Qgf 
for the lower mass bin. Overall, it is remarkable that regardless of 
methodology or feedback processes, all models predict a^i to be 
around unity with little evolution from z = 4 ^ 0, regardless of 
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Figure 2. Top panel: Star formation activity pai'ameter a^f = 
(A/*/SFR)/(tH - IGyr) at ^ = - 4, evaluated at M, = lO^O-'^Afo. 
Red circles show results from hydro simulations with a momentum-driven 
outflow model, cyan circles {z ^ 2) show hydro simulations with no out- 
flows, and magenta circ les (2^2 show "constant wind" outflow simu- 
lations (see Oppenhei mer & Davell2006i . for details of outflow models). 
Green triangles sho w the Millen i um SA M results, and blue crosses show 
the observation s bvlElbaz et al. I i2007l z ^ 0, 1), iNoeske et afl j2007al . 
2 Si 0.45), and lDaddi et al. I )2007l 2 S 2). In aU cases, regardless of mod- 
eling technique or outflow implementation, cigf is predicted to be around 
unity or larger at all epochs and at all masses. The observed evolution of 
Ogf out to 2 ^ 2 is markedly discrepant from model predictions. Bottom 
panel: Same as above for = 10'^'' Mq, showing that the trends are 
broadly insensitive to stellar mass. 



stellar mass. This suggests that a^i ~ 1 is a generic consequence of 
smooth and steady cold accretion dominating the growth of these 
galaxies. 

In contrast, the observed Af, — SFR relation yields a rapidly 
rising a^i = 0.19 0.31 -» 0.75 from 2 = 2-»l-^0at 
IO^^ ^'A/q. Values at lO^'Af© are slightly lower but show simi- 
lar rapid evolution. The dramatic difference between observed and 
predicted amplitude evolution either indicates a significant misun- 
derstanding in our basic picture of galaxy assembly, or else that 
observations have been misinterpreted in some way. Given the cur- 
rently uncertain state of galaxy formation theory, one might be in- 
clined towards the former, and call into question the entire picture 
of cold mode-dominated accretion. However, it turns out that the 
rapidly declining a^f also results in uncomfortable conflicts with 
other observations of high-redshift galaxies. To see this, let us dis- 
card the idea of cold mode accretion, or any preconceived theoreti- 
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cal notion, and ask: What possible star formation histories for these 
galaxies could yield Osf rising with time as observed? 

One scenario is to postulate that observed high-z galaxies are 
progressively more dominated by bursts of star formation. Indeed, 
some models have suggested that bursts may c omprise a significan t 
portion of the hig h-z galaxy popu lation (e.g. iKolatt et al. Ill99si) . 
However, as Noesk e et al. I j2007ah points out, this seems unlikely 
given the tightness of the Af,-SFR relation. They argue that the 
observed scatter of 0.3 dex at z ~ 1 implies that most star form- 
ing galaxies have current star formation rates < x 2 — 3 of their 
quiescent value, and that the mean of the relation represents the 
quiescent "main sequence" of galaxy assembly. Few galaxies be- 
tween this main sequence and the passive population are seen, ar- 
guing for rapid movement between the t wo populations. At z ~ 2, 
the constraints are even more severe, as iDaddi et al. I ilOOH) mea- 
sure a la scatter of 0.29 dex (0.16 dex interquartile), while the 
observed asi ~ 0.2 implies they are bursting at ~ x5 their quies- 
cent rate. Since these observations are able to probe SFRs down to 
~ 10 Af0/yr, moderate stellar mass (M. > 10^°-^ Mq) galaxies 
with Qfsf ~ 1 would have been seen if present, but are not. The 
only way the burst scenario would work is if all observed galaxies 
somehow conspire to be bursting at approximately the same level, 
which seems highly contrived. Hence bursts do not seem to be a 
viable way to produce low a^f . 

Another possible scenario is that models simply need to be- 
gin forming stars at a later epoch. This would effectively lower 
Qaf as defined here, while still maintaining a fairly constant SFR 
as required by the observed scatter. Theoretically, such a sce- 
nario is easier said than done, as the epoch where galaxy forma- 
tion begins is not a free parameter even in semi-analytic mod- 
els, but rather is governed by halo formation times that are well- 
constrained by current cosmological parameters. Such a scenario 
would require overwhelmingly strong feedback at early epochs, 
far stronger than even the strong feedback present in e.g. hydro 
simulations t hat match ob served z ~ 4 — 7 ga laxy star forma- 
tion rates jPave. Finlator & Oppenheimerl l2006l : iSouwens et al. I 
l2007h , in order for copious dense gas to be present in early halos 
yet not form stars. Still, in the spirit of discarding any theoretical 
preconceptions, "delayed galaxy formation" appears to be a viable 
solution. 

However, quantitative consideration of delayed galaxy forma- 
tion leads to significant problems when confronted by other obser- 
vations. Al z — 2, Qfsf ~ 0.2 implies that (constant) star formation 
began roughly 0.5 Gyr prior ~ which corresponds to jz ~ 2.3! Of 
course it cannot be true that all galaxies began forming stars at 2 ~ 
2.3, as many star-forming galaxies are seen at significantly higher 
redshifts (e.g. 'Bouwens et al. "2007"). Indeed, both the cosmic star 
formation rate (e.g. Hopkins & Beacom 2006; Fardal et al. 2007), 
and the rest-UV luminosity function (Reddv et al. II2007I) are ob- 
served to be fairly constant between z ~ 4 ^ 2. So it appears 
galaxies were forming stars quite vigorously at z ^ 2.3. Even 
if galaxies smoothly ramped up their star formation to values ob- 
served at 2 ~ 2, they would still need to start at z ~ 2.7. 

If most galaxies cannot have such a low Ogf, perhaps the ob- 
served systems represent only a small subset of galaxies at 2: ^ 2, 
and the rest are simply unseen. Since current surveys identify both 
UV-luminous and red, dusty star forming galaxies, it must be that 
the unseen galaxies are passive, or else so extremely extincted that 
they only show up in e.g . sub-millimeter surveys. For the latter 
case, iReddv et al. [ fcOOSi) found that the cosmic SFR has a rela- 
tively minor contribution from such severely extincted systems that 
escape BzK/LBG selection at 2 ~ 2, so they are not likely to be the 



dominant growth mode for galaxies. Moreover, this actually exac- 
erbates the M, —SFR problem, because it provides a mode to form 
many more stars beyond what is formed quiescently, which would 
make the expected stellar masses even higher than just from the 
quiescent mode. 

Hence the unseen galaxies must be passive, or forming stars 
at some und etectably low rate. In the currently popular downsiz- 
ing scenario ('C owieetal.ll9 96l). these systems would have formed 
their stars earlier in a rapid mode, and have now attenuated their star 
formation through some feedback process so as not to be selected 
in star-forming samples. Perhaps by 2 ~ 2, downsizing is already 
well underway, and the observed galaxies only represent the small 
fraction of systems that have just recently begun forming stars and 
will soon be quenched. 

While qualitatively sensible, the "unseen passive galaxies" 
model is once again quantit atively inconsistent wi th other obser- 
vations of high-2 galaxies. IPerez-Gonzalez et al. ] |2007) used a 
5p//zer/IRAC-selected sample to determine that the star-forming 
BzK (sBzK ) sele ction technique, i.e. the one employed by 
iDaddi et al. I i2007h to measure M«-SFR at 2 ~ 2, encompasses 
88% of their galaxies (2691/3044) from 1.6 < 2 < 2.5. This is re- 
markably efficient, significantly more so than e.g. the Lyman break 
or Distant Red Galaxy selection techniques in obtaining a stellar 
mass-limited sample. Their stellar mass completeness is estimated 
to be ~ 75% above M* > lO^^'^M© for passively evolving galax- 
ies. Putting these together, this means that conservatively at least 
two-thirds of all galaxies above this mass threshold are selected 
as sBzK galaxies, and hence less than one-third are unseen pas- 
sive galaxies. In contrast, the observed M,— SFR relation at 2 ~ 2 
implies Qsf ~ 0.2, meaning that ~ 80% of galaxies should be 
unseen (passive), or else the duty cycle of such objects should be 
~ 20%. This assumes galaxies began significant activity ~ 1 Gyr 
after the Big Bang, and that Qsi — 0.2 at 2 > 2; the latter is again 
conservative, as an extrapolation would indicate that Ogf should be 
even lower at z > 2. Hence it does not appear that sBzK selection 
misses nearly enough systems to make the unseen passive galaxies 
scenario viable to explain Ogf (2 ~ 2). 

Another way to accomodate the low scatter, unity slope, and 
rapid evolution in A/, —SFR is to propose that ga laxies evolv e 
along the A/*— SFR relation, as proposed by Daddi e t al. I ( |2007|) . 
At 2 ~ 2, Osf = 0.2 implies that galaxies must be growing very 
quickly. Rapid growth means that for any given galaxy, the am- 
plitude of Af*— SFR is essentially unchanging. Now since SFRoc 
A/, = J SFR dt, an unchanging amplitude implies a solution of 

SFR = Cexp(f/r) =^ A/, = Cr exp(t/r) = r SFR. (2) 

In other words, in this scenario, galaxies have exponentially grow- 
ing star formation rates and stellar masses ! This growth would pre- 
sumably continue until star formation is truncated by some feed- 
back mechanism, probably when Al, > 10^^ Af©. While at first 
this "exponentially growing galaxies" idea may seem incredible, 
of all the scenarios it actually comes closest to reconciling all the 
high-2 galaxy data. 

One way to test this scenario is to consider the evolution of the 
number densities n of galaxies. The e-folding time of stellar mass 
growth is given by r = A/, /SFR. Assuming (conservatively) that 
Qsf ~ 0.2 remains constant at 2 > 2, then r oc t„ - 1 Gyr0 



^ Technically, the fact that r changes with time means that the assumption 
of a constant Af, -SFR amplitude in equation|2]is invalid; however, a simple 
evolutionary model accounting for this yields virtually identical results. 
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A simple calculation shows that there are e.g. five e-folding times 
(or 2.2 decades of growth) from z = 3.5 2. For a given stellar 
mass M*, n(M,) at z; = 2 should then b e similar to n(M ,/10^-^) 
at z = 3.5 (neglecting mergers). Font ana et al. 1 ( l2006l) measured 
the number density of galaxies with M, = 10^^ at 2: ~ 2 to 
be 7 X 10~*Mpc~^, and those with M. = IO^ ^Mq a.t z = 3.5 
have n — 3 X 10^'^Mpc^'^ (based on an extrapola t ion of their 
STY fits), and stellar mass functions from lElsner et al. yield 
similar values. Hence at face value current data suggests that the 
stellar mass function does not evolve quite as rapidly as expected in 
this scenario. Still there are substantial uncertainties in the observed 
stellar mass function at such low masses and high redshifts, so a 
factor of a few discrepancy is not compelling enough to rule out 
the exponentially growing galaxies scenario. 

Another scenario that has been shown to repro duce M^-SFR 
evolutio n is the staged galaxy formation model of 'Noeske et al.J 
( l2007bl) . Here, a delay is introduced in the onset of galaxy forma- 
tion that is inversely related to the galaxy mass. Hence large-mass 
objects form early and quickly, while small ones start forming later 
over longer timescales. This is similar to the exponentially-growing 
galaxies scenario in the sense that the latter implicitly assumes a 
sort of delay, because galaxies must accumulate lots of gas in their 
halo without forming stars in order to have a sufficient gas reser- 
voir to then consume exponentially. Hence it may be possible to 
constrain staged galaxy formation with similar measurements of 
the stellar mass function to higher-z; this is being investigated by 
K. Noeske (priv. comm.). If either exponentially-growing or staged 
galaxy formation is correct, this would indicate that star forma- 
tion at early epochs d oes not scale with gas density as seen lo- 
cally jKennicuttll"99^ , or else that cold gas is somehow prevented 
from accumulating in halos as expected from growth of structure. 

In summary, the small Ogf at z ~ 2 implied by the A/, — SFR 
amplitude is difficult to reconcile with the observed scatter in 
M*— SFR, the observed evolution of star-forming galaxies to 
higher redshifts, and the observed fraction of passive galaxies. 
Taken together, the latter data all point to the idea that observed 
star-forming galaxies represent the majority of galaxies at interme- 
diate masses 10^" — lO^^A/0) at 2 ~ 2, and have been qui- 
escently forming stars over much of a Hubble time, i.e. a^f ~ 1. 
While designer models of stellar assembly cannot be excluded, they 
would be highly unexpected, and would require substantial revision 
to our understanding of how gas accumulates and forms stars in hi- 
erarchical structure formation. The general difficulty in reconciling 
the data together with the models' preference for asf ~ 1 suggests 
an uncomfortable situation that bears investigation into alternative 
solutions. 



4 SYSTEMATIC EFFECTS IN AL AND SFR 
DETERMINATIONS 

Determining Al, and SFR from broad-band observations is fraught 
with a notoriously large number of systematic uncertainties that 
could alter the values of a^f inferred from observations. Here 
it is considered whether such systematics could be responsible 
for the low value o f Osi at high-«, particularly focusing on the 
iDaddi et al. I ilOOt analysis at z ~ 2 since that is where the 
discrepancies from model expectations and the range of possi- 
ble systematics are greatest. Star formation rates there are gener- 
ally estimated by summing rest-UV flux (uncorrected for extinc- 
tion) with 24fi flux. For objects that are anomalously high or low 
in 24/i to dust-corrected UV SFR estimations, they use the dust- 



corrected UV or 24/i estimate, respectively; however, only one- 
third of their sample are anomalous in that sense. Meanwhile, stel- 
lar masses are estimated from SED fitting: Using multi-band pho- 
tometry or spectra, one fits model s tellar population templates (in 
their case, iBruzual & Chariot |2003l). varying par ameters until the 
best fit is obtained. In detail. iDaddi et al. I ( |2007|) used a calibra- 
tion ofBzKcolortostellar mass based on fitting K20 survey spec- 
tra ( Fonta na et al. 1120041). The key assumptions in SED fitting are 
population synthesis models, dust extinction models, and the form 
of the star formation history, all of which introduce systematic un- 
certainties. 

Current uncertainties in obtaining accurate stellar masses from 
rest near-IR data revolve around the poorly known contribution of 
thermally-puls ating asymptotic giant branch (TP-AGB) stars to the 
if-band light jMaraston et al. II2OO6I) . This is found to lower the 
stellar mass estimates for high- 2: galaxies by a factor of two or 
more (Maraston et al. 2006; Bruzual 2007). However, the sense 
would be to exacerbate the problem b y lowering Ad, and hence 
further reducing as t, compared wi t h the l Bruzual & Chariot 1l l2003h 
models utilized in iFontana et afl {20o57 Hence while uncertain- 
ties in population synthesis models are not negligible, they seem 
unlikely to make a > x 3 difference in the proper direction. 

It is possible to hide more stellar mass in objects if 
one postulates a very old underlying stellar population, 
in other words by modifyin g the assumed SFH form. 
IPapovich. Dickinson & Ferguson I ( 1200 ih found that for UV- 
selected galaxies, a mass increase of ~ x3 — 8 was possible 
without violating r est-optical constrai nts. However, 7f-selected 
samples like that of iDaddi et al. I ^20071) are generally redder with 
more contribution from old stellar light, so the amount that one can 
hide is more limited (C. Papovich, priv. comm.). More problematic 
is the fact that doing so requires assuming that the bulk of the 
galaxy's stars formed very early in the Universe, and hence the 
current SFR is much lower than at that early epoch. It turns out 
that this is not a self-consistent solution for reducing ct^i- This is 
because if such a SFH history were true, one would expect that the 
current M* would be significantly more than obtained by assuming 
a constant SFH over a Hubble time. This would imply a^t S> 1 
for those galaxies. Hence even hiding 5x as much stellar mass 
as currently inferred would be insufficient, because it would only 
make ^ 1, whereas such a SFH implies a^i 3> 1. In order 
for there to be a large old stellar population while still having a 
high current SFR, it must be that such galaxies formed many stars 
early, then sat around not forming stars until a recent flare-up; 
however, this scenario suffers the same difficulties as the delayed 
galaxy formation model. Hence one cannot self-consistently alter 
the SFHs to hide a large amount of very old stars in these galaxies. 

Instead, it may be possible to appeal to some systematic re- 
duction in the inferred star formation rate. Why might the true SFR 
be lower than estimated? One possibility is that the extinction cor- 
rections being used for UV-estimated SF Rs are too excessive . Now 
recall that for the majority of galaxies, iDaddi et al. I ( l2007 f) esti- 
mated the SFR from the uncorrected UV plus 24/i flux, so for those 
cases the extinction correction is not relevant; it is only relevant 
for the < 25% that are "mid-IR excess" systems. For those, there 
may be room to alter the extinction law such that t he extinction 
inferred from the UV spectral slope is overestimated. iDaddi et al. I 
[2007), like most stu dies of actively star forming galaxies, use the 
ICalzetti et aTl ( |2000|) law. Using instead a Milky Way or LMC law 
would indeed yield less extinction for a given UV spectral slope, 
but the typical reduction is less than a factor of two. So it seems 
unlikely that extinction alone is causing the low Qsf even in the mi- 
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nority mid-IR excess objects, unless the extinction law were dra- 
matically different than anything seen in the local universe. 

Another possibility is that there is an additional contribution 
to the rest-UV flux besides star formation. The obvious candidate 
are AGN, which are co mmon in 2: ~ 2 star forming galaxies. While 
IPapovich et Z1 j200d ) found that only ~ 25 % of > 10"A^a 
galaxies at 2 ~ 2 have X-ray detectable AGN. IOaddi et al. iJioOTbh 
identified obscured AGN through their mid-IR excess, and found 
that they appear in ~ 50% of galaxies at A/* > 4 x 10^° M©, and 
dominate the overall AGN population at z ~ 2. I mportantly, they 
also found (in agreement with previous studies, e.g.^ Papovich et alH 
[2OO6) that smaller m ass galaxies have a lo wer AGN fr action, as ex- 
pecte d theoretically faopkins et al. |2007|) . Recall that lDaddi et al. I 
( 120071) used UV fluxes to estimate SFR in such mid-IR excess ob- 
jects, and by comparison with X-ray data they estimated that such 
an inferred SFR is unlikely to be wrong by a significant factor. 
Moreover, since AGN are increasingly rare in smaller mass sys- 
tems, it would be an odd coincidence if they mimicked the relation 
M« ocSFR'''^ so tightly. Like the burst model, this represents a 
fine-tuning problem, where AGN have to contribute a very specific 
amount to the inferred SFR at each mass scale for no obvious phys- 
ical reason. This cannot be ruled out, but seems a priori unlikely. 

It is also possible that AGN could be present in galaxies that 
are not mid-IR excess systems. The strong dichotomy in the X-ray 
properties suggest there is not a continuum of AGN, and that AGN s 
are limited to mid-IR excess systems (e.g. Daddi et ^TlbOOTbh . 
However, there exist AGN templates where the contibution to rest- 
UV and r est-8/.t flux are br oadly comparable to that in star-forming 
galaxies jPale et al. Il200d) . These could be mistaken for "normal" 
star-forming systems, whereas in fact the dominant flux is from 
AGN. Still, once again, to mimic the tight —SFR relation would 
require AGN to be contaminating the SFR estimate at a similar 
level at all stellar masses. Moreover, if they are dominating the UV 
flux then they should appear as central point sources, but Hub- 
ble/ACS imaging of sBzK galaxies generally do not show such 
bright central sources (M. Dickinson, priv. comm.). Hence while 
AGN are certainly present, it seems unlikely that they are a major 
contaminant for UV-derived SFRs. 

A perhaps even more troublesome issue is the use of rest-8/i 
(i.e. 24/1 at 2 ~ 2) as a star formation rate indicator. This is quite 
uncertain, because the dominant flux at those wavelengths is from 
polycyclic aromatic hydrocarbon ( R'VH) emission feat ures, whose 
nature is not well understood. As ISmith et al. I ilOOH) points out, 
there is more than a factor of two scatter in star formation rate indi- 
cators that assume a fixed (local) PAH template. Moreover, owing 
to PAH features moving in and out of bands, standard photometric 
redshift uncertainties of Az/(1 + 2) ~ 0.1 can result in up to an 
order of magnitude error in inferred SFRs. On the other hand, the 
tight scatter in the 2 ~ 2 Af»— SFR relation would seem to limit 
the amount of scatter owing to PAH feature variations, or else be an 
indication that UV light dominat es over PAH emiss i on (though the 
latter is not what is seen locallv: ISmith et aL 1l2007l) . ICaizettietal.l 
( I2OO7I) notes that 8/i emission has a substantial metallicity depen- 
dence, but it is in the sense that lower metallicities have lower emis- 
sion, so their SFR would be underestimated by using solar metal- 
licity calibrations. So while it is conceivable that some unknown 
reason causes PAH emission per unit star formation to be system- 
atically higher in high-2 galaxies compared to local ones, thereby 
causing the SFR to be overestimated, but there are no obvious indi- 
cations that this would mimic a tight Af»-SFR relation. 

Finally, the observed SFR and A/, 's in lDaddi et al. I ( 2007h are 
inferred assuming a Salpeter IMF. Assuming instead a Chabrierl 



( I2OO3I) or lKroupalfcoOlh IMF would result i n shifting both AL and 
SFR down by « 0. 15 dex felbaz et al. l2007h . Hence to first order it 
makes no difference as long as the Af»— SFR slope is around unity 
(E. Daddi, priv. comm.). Fundamentally, this is because at high- 
2 both stellar mass and star formation rate indicators are driven 
by light from stars > 0.5A/q, and favored present-day Galactic 
IMFs generally agree on the shape above this mass. Hence more 
substantial IMF variations would be required. 

In summary, there are no obvious paths to systematically alter 
Aft and SFR determination to make up up the required ~ x 3 — 5 
difference in a^f at 2 ~ 2. It is true that many of these systematics 
cannot be ruled out entirely, and it may be possible that a combina- 
tion of such effects could explain part or even all of the difference. 
It is also possible that locally calibrated relations to estimate Af* 
and SFR may be substantially different at high-2 for unknown and 
unanticipated reasons. The reader is left to assess the plausibility of 
such scenarios. The view taken here is that the difficulty in obtain- 
ing a straightforward solution motivates the consideration of more 
exotic possibilities. 



5 AN EVOLVING IMF? 

One possible way to alter A/, -SFR evolution is to invoke an IMF 
that is in some direct or indirect way redshift-dependent. The key 
point is that most direct measures of the star formation rate ac- 
tually trace high-mass star formation, while stellar mass is domi- 
nated by lower-mass stars. Hence by modifying the ratio of high 
to low mass stars formed, i.e. the shape of the IMF, it is possible 
to alter the A/, —SFR relation. The possibility of a varying IMF 
has been broached many times in various contexts (e.g. lUarsonl 
ll998l ; lFerguson et al. l20ol : lLacev et al. l2007l;lFardal et al. l2007hT 
though typically as a last resort scenario. It is highly speculative as 
there is no clear-cut evidence at present that supports a time- or 
space- varying IMF, but it is worth investigating in light of growing 
observa tional and theore tical arguments in its favor. 

As iKroupa 1 l l200lh points out, a universal IMF is not to be 
expected theoretically, though no variations have been unequivo- 
cally detected in local studies of star-forming regions. The aver- 
age Galactic IMF appears to be well-represented by several power 
laws, steepest at > lAf©, then significantly shallower at < 0.5Af(7), 
with a turnover at < O.lAf© (see review by Kroupa ^l 2007l h. Forms 
that follow this include the currently-favored Kroupa J ( 120011) and 
IChabrierl (2003) IMFs. But deviations are seen, at least at face 
value. The nearby starburst galaxy M82 appears to have an IMF 
deficient in low-mass stars dRieke et al. Ill98(i Il993l) . A similarly 
bottom-light IMF is inferred in the highly active Arches cluster near 
the Galactic center, sug gesting that vigorous star formation makes 
the IMF top-heavy (e.g. iFigerl2005h . Also, the youn g LMC cluster 
R136 shows a flattening in the IMF below ^ 2Mq { Sirianni et al. I 
I2OOOI) . There are hints that presently active star clusters form 
more low-mass stars compared to the disk-averaged 5 Gyr-old 
IMF (Ki'oupa 2007), suggesting that the IMF was weighted to- 
wards heavier stars at lower metallicities and/or earlier times. Even 
if the IMF shape is universal, larger starforming regions may pro- 
duce more high-mass stars simpl y because they have enough mate- 
rial to aggregate into larger stars dWeidner & Kroupa Il200d) . 

Hence there are hints of IMF variations, albeit controversial, 
and interestingly they consistently go in the direction of having 
a higher ratio of massive to low-mass stars in conditions sim- 
ilar to those in high-2 galaxies. Compared to present-day star- 
forming systems, galaxies at high-2 have far higher star formation 
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surface densities than the Galactic disk l lErbetal. 1120063) , have 
higher gas content and hence presumably more massive starform- 
ing clouds (Srb et al. 2006b), and have somewhat lower metallic- 
ities dSavaglio et al. 2005; Erb et al. 2006a). Larson (1998) and 
iFardal et al. ( 2007 ) list other pieces of circumstantial evidence that 
favor early top-heavy IMFs. So while there are no direct observa- 
tional evidences or strong theoretical arguments for it, it is not in- 
conceivable that the IMF in high-2 galaxies may be more top-heavy 
or bottom-lighjfl than the present-day Galactic IMF. 

Such an IMF would have three main effects on the M*— SFR 
relation: It would increase the output of UV flux per unit stellar 
mass formed; it would cause more stellar mass loss from stellar 
evolution; and the mass recycling would provide a larger gas reser- 
voir with which to form stars. All three effects go towards recon- 
ciling theoretical expectations with observations of the i\f*-SFR 
relation. The SFR inferred by assuming a standard IMF would be 
an overestimate, since many fewer low-mass stars are formed per 
high-mass star. Increased recycling losses would lower the theoret- 
ical expectations for the amount of stellar mass remaining in these 
galaxies, bringing them more into line with observed stellar masses. 
Recycled gas would increase the available reservoir for new stars, 
effectively delaying star formation in galaxies and pushing models 
towards the "delayed galaxy formation" scenario. Hence even mod- 
est evolution in the IMF could in principle reconcile the observed 
and theoretically-expected a^i evolution. 



6 A SIMPLE IMF EVOLUTION MODEL 

How is the IMF expected to evolve with time? This is difficult to 
predict since no a b initio theory for the I MF exists today. But a sim- 
ple conjecture bv iLarson] jl998ll200^ may have some empirical 
value. He notes that there is a mas s scale of ^ O.SMq below which 
the Galact ic IMF becomes flatter ( K roupa, Tout. & Gilmore Il993l ; 
iKroupa 1200 1 ; Chabrier 2003), and at which point the mass contri- 
bution per logarithmic bin (MdN/d log M) is maximized. Larson 
uses simple thermal physics arguments to show that this charac- 
teristic IMF mass scale (A/imf) is related to the minimum tem- 
perature Tmin of molecular clouds, measured today to be around 
8 K. In that case, if the temperature in the interstellar media (ISM) 
of high- 2; galaxies is higher, this characteristic mass will shift up- 
wards, resulting in an IMF that is more bottom-light. 

To generate a crude IMF evolution model, let us assume that 
the i\fiMF oc (1 + z)^, where 7 is some u nknown param eter whose 
value is set by the various effects above. iLarson 11 1 19851) predicted 
analytica lly that Mimf oc rj ,;^^, while recent numerical experi- 
ments bv lJapDsen et al. I ( |2005|) found a less steep scaling of T^jn- 
So for example an evolution that scaled purely with CMB tem- 
perature would yield 7 = 1.7. However, other factors may cause 
an increased ISM temperature, e.g., elevated disk SFR rates caus- 
ing more supernova heat input; lower cooling rates owing to lower 
metallicities; and/or photoionization from a stronger metagalactic 
UV flux (e.g. lScottet al. 2002). 

The approach adopted here is to treat 7 as a free parameter, 
to be constrained using the observed evolution of a^i. The assump- 
tion is that asi{z — 0) reflects the present-day Galactic IMF, while 



^ The distinction between top-heavy and bottom-light IMFs is mostly se- 
mantic. Here the convention is used that top-heavy means a high-mass slope 
that is less steep than the local Salpeter value, while bottom-light means that 
the high-mass slope remains similar but low-mass stars are suppressed. 



asi{z > 0) is inferred to be lower purely because of the assump- 
tion of a non-evolving IMF, whereas in reality Osf («) (at any given 
mass) is constant as expected in models. 

An "evolving Kroupa" IMF is defined in which 



Mimf = 0.5(1 + ^)^M0. 



(3) 



Above Mimf, the IMF has a form dN/dlogM cx Af"^'^, while 



below it scales as M 



. Note that the current Galactic IMFs also 



have a turn-down at masses below 0.1 Mq ; the present analysis pro- 
vides no constraints in this regime, as the stellar mass formed there 
is a small fr action of t he tot al. The exact form of the IMF is not 
critical here; I Kroupa 1 j200lh is chosen for its convenient explicit 
parameterization in terms of i\/iMF, but the analysis below could 
equally well have been done with the similar Chabrier IMF. 

The goal of the evolving IMF is to increase the ratio of the star 
formation rate of high-mass stars (SFRhim) to total stellar mass 
formed, such that the inferred Qsf(z) is non-evolving. Note that 
different probes such as UV luminosity, Ha emission, and mid-IR 
emission trace slightly different regimes of high-mass star forma- 
tion, but so long as the high-mass IMF slope remains fixed, the 
relative ratios of such emission should be constant. So for exam- 
ple, in order to produce no evolution in a^i from 2; ~ —* 2 at 
IO^^Mq, the factors by which / =SFRhiM/A^» must be raised are 
1.9 at 2 = 1 and 3.3 at z = 2, relative to a standard Kroupa IMF 
(cf. Figure[2]l. Hence 7 must be chosen to match these factors. 

What value of 7 wo uld produce such evolution? To deter- 
mine this, the PEGASE.2 iFioc & Rocca-Volmerangelll997i) pop- 
ulation synthesis model is employed, taking advantage of its fea- 
ture allowing user-definable IMFs. PEGASE outputs a variety of 
quantities as a function of time from the onset of star formation, 
given an input star formation history. The assumed SFH is taken 
to be exponentially-declining with a 10 Gyr decay time and start- 
ing 1 Gyr after the Big Bang, in order to approximate a typical 
model SFH. An example of such a SFH is shown as the magenta 
line in the lower right panel Figure[T] normalized to a star formation 
rate of 2 Afoyr^^ today; as is evident, it is quite similar to typical 
SFHs from the simulations, at least for moderate stellar masses. A 
constant SFH was also tried, and the results were not significantly 
different. Besides the input SFH, all other parameters are taken as 
the PEGASE default values. No dust extinction is assumed. 

The PEGASE output quantities specifically used here are the 
stellar mass Af, formed, and the un-extincted 1650A luminosity 
(ii65o)- The latter is used as a proxy for SFRhiM. Using these, it 
is possible to compute / for any given assumed IMF at any given 
time. Note that PEGASE accounts for both stellar mass accumula- 
tion and mass loss owing to stellar evolution, but does not account 
for the impact of recycled gas. By varying A^imf, it is possible to 
estimate / for any assumed 7, and thereby determine what value of 
7 would yield no a^f evolution. 

In practice, determining / =SFRhiM/M« at any epoch re- 
quires integrating the mass growth up to that epoch for an evolving 
IMF. To do so, first PEGASE models are computed for a wide range 
of A/iMF values. Then the following procedure is used: 

(1) A value of 7 is selected; 

(2) the following integral is computed: 



M,(to)=/ -^{MjMF,t)dt, 



(4) 



where AM / At is the stellar mass growth rate at time t taken from 
PEGASE models interpolated to the appropriate A/imf as given in 
equation[3l 

(3) 1/1650 is taken at time t from PEGASE, and /ov(to) = 
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Figure 3. Evolution of the ratio of high-mass star formation rate to stel- 
lar mass, for the evolving Kroupa IMF (equation [5) relative to a standard 
Kroupa IMF, normalized to unity at 2 = 0. The curves are computed from 
PEGASE.2 models for different values of 7, as described in ^ The data 
points are from fits to the observed A/» — SFR relation at 2 = 0.45, 1, 2 (de- 
scribed and shown in Figure|2j at stellar masses of IO^^'^Mq (squares), 
10^*^ Mq (triangles), and 10^'^ Mq (circles). The required evolution is 
reasonably well fit by 7 = 2, i.e. Mimf = 0.5(1 + zf^Q. Long 
and short dashed-dot models show results for the iLarson I jl998l) scenario 
where the minimum temperature of molecular clouds is set by the C MB 
temperatur e, and Mimf 

t". , where /3 = 3.35 jLarson |[T983) and 
/3 = 1 . 7 I JaDDsen et al. 2005), respectively. Because the CMB temperature 
does not exceed the the local minimum temperature in molecular clouds of 
8 K until 2 > 2, this scenario does not yield substantial evolution at 2 < 2. 



ii65o(io)/A^, (to) is computed (the subscript "ev" stands for 
"evolving"); 

(4) Similarly, fst{to) is computed for a standard Kroupa IMF, 
where A/imf = O.SM©; 

(5) The ratio /ev(to)//st(to) is computed, and compared to the 
factors needed to produce no asi evolution. 

The results of this procedure are shown in Figure [S] namely 
/ev//st as a function of redshift for several values of 7. The val- 
ues by which the A/, —SFR relation needs to be adjusted to pro- 
duce no ftsf evolution at z — 0.45, 1, 2 are indicated by the sym- 
bols: Squares are for M* = 10^'^ Mq, triangles for lO^^M©, and 
circles for lO^^'^Afp). These are computed from published fits to 
M«- SFR i lNoeske etal.1l2007al : lElbaz et al. I bool iDaddi et al. I 
l2007h . Because the Af»— SFR relation becomes shallower with 
time, the required adjustment factors are generally smaller at lower 
masses. 

Figure [3] shows that 7 = 2 provides a reasonable fit to the 
required IMF evolution for typical star-forming galaxies at A/, — 
10^°Mq. The impUed values of A^imf are 2A'/q and 4.5M0 at 
z = 1,2, respectively. Note that there is no a priori reason why a 
single value of 7 should fit at all redshifts considered, so the good 
fit supports the form of A/imf evolution assumed in equation [3] 

As a check, an independent method for determining A/imf 
evolution was explored that is less comprehensive, but conceptu- 
ally more straightforward. Here, the high-mass star formation rate 
is assumed to be unchanged, while observations of stellar mass 
are assumed to reflect stars near the main sequence tumoff mass 
(A/turnoff), since red optical light is dominated by giant stars. For 
an evolving Kroupa IMF, to obtain a reduction of a factor of x in 
the amount of stars produced at Mtumoff, it is easy to show that 
log AfiMF = log Aftuinofi + log X. A/turnoff may be estimated by 



noting that stellar lifetimes scale as M , and that the Sun has a 
lifetime of 10 Gyr; this yields Mt^.^^ff = 10/(tH - 1), (tn in 
Gyr), which is 1.24M0 and I.SSA/q &i z = I, 2, respectively. In- 
serting values of a; = 1.9,3.3 at z = 1,2 yields A/imf = '2.A 
and 5.2AfQ. These values are quite similar to those obtained from 
PEGASE modeling, showing that the results are not critically de- 
pendent on details of PEGASE. 

The same pr ocedure can be applied to the scenario proposed 
by iLarson I ( Il99^) . where the CMB temperature is solely respon- 
sible for setting a floor to the ISM temperatur e. In that case , 
Tmiu = MAXfTcMB.SK] and Mimf o c , . [larsonJ <tl98^ 
predicted /3 = 3.35, while Ijappsen et al. I liooa) found /? = 1.7. 
The results of this scenario for these two values of [3 are shown as 
the long and short dot-dashed lines, respectively, in Figure|3] As is 
evident, this IMF evolution is not nearly sufficient to reconcile the 
observed and predicted Usi evolution out to z ~ 2, mainly because 
TcMB only exceeds 8 K at 2; > 1.93, and hence A/imf does not 
change from 2; ~ — 2. 

Note that since CMB heating of the ISM is subdominant at 
all redshifts, the assumed form of IMF evolution (eq. [3] scaling as 
^cmb) is rio' physically well motivated. While the evolution out to 
z ~ 2 is well-fit by such a form, it could be that the form changes at 
higher z, or that it should be parameterized as some other function 
of z. For instance, if it is the vigorousness of star formation activity 
that determines A/imf, perhaps IMF evolution actually reverses at 
very high redshifts. For lack of better constraints, the form for IMF 
evolution in equation [3] will be assumed at all z, but this should be 
taken as an illustrative example rather than a well-motivated model. 

This analysis also makes predictions for Ofsf at 2; > 2 that 
would be inferred assuming a standard (non-evolving) IMF. For in- 
stance. Figure |3] implies that at IO^^A/q, asf (z = 3) = 0.13 and 
asi{z = 4) = 0.12. Notably, the evolution slows significantly at 
z > 3, and actually reverses at 2; > 4, though the simple SFH as- 
sumed in the PEGASE modeling may be insufficiently realistic to 
yield valid results at 2; > 4. Any such extrapolation of this IMF evo- 
lution should be made with caution, as the A/* -SFR relation only 
constrains it out to z ~ 2. However, in the next section this will 
be compared to observations that suggest that such IMF evolution 
may be reasonable out to 2; ~ 3+. 

Figure |4] illustrates the evolving IMF, in terms of stellar mass 
formed per unit logarithmic mass bin. A Salpeter IMF is shown 
for reference, and evolving Kroupa IMFs are shown at z = 0, 1, 2 
with characteristic masses A/imf ~ 0.5Mq, 2Mq, and 4.5A/0, 
respectively. As expected, the evolution results in more high mass 
stars compared to low mass ones at higher redshifts. While such a 
dramatic increase in A/imf may seem surprising, recall that stud- 
ies of local highly active star-fo rming regions suggest a character- 
istic mass scale of a few Mq dRieke et al. 1 1 19931 ; ISiriaimi et al. I 



l200(]| ; lFigerll2005l). Note that an IMF with a steeper high-mass 



slope (e.g. lMiller & Scalolll979l : lKroupa. Tout. & Gilmorelll993h 



wo uld result in eve n stronger IMF evolution being required than in 
the lKroupal ( 1200 ih case. 

In summary, an evolving Kroupa IMF whose characteristic 
mass increases with redshift as A/imf = 0.5(1 + z)^Mq is able 
to alter the observed a^i evolution from 2; ~ 2 — ^ into one 
with approximately no evolution, as predicted by models or in- 
ferred from complementary observations. The precise IMF shape 
is not well constrained; it is certainly possible that an IMF with 
different behavior, for instan ce an evolving high-mass slope (e.g. 
iBaldrv & Glazebrookll2003h . could produce the same Qgf evolu- 
tion. The only requirements are that it produces the correct ratio of 
high-mass stars relative to low-mass stars as a function of redshift. 
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Figure 4. Initial mass functions, multiplied by M to show the mass fomied 
per unit loga rithmic m a ss bi n, and arbitrarily normalized at 1 Mq . Salpeter 
(dashed) and lKroupal EoOl) ( green solid) IMFs are shown; Chabrier is 
similar to Kroupa. Evolving Kroupa IMFs are shown at z = 1 (magenta 
got-dashed) and z = 2 (red dashed), with characteristic masses Mimf = 
2Mq and 4.5Mq, respectively, as required to remove evolution in a^t 
(see ij6). The cumulative z = IMF (cyan dotted) shows the summation 
of an evolving Kroupa IMF with Mimf = 0.5(1 + z)'^ over the history 
of the Universe, assuming an evolving IMF cosmic star formation history 
as shown in Figure [5] Also shown is the "paunchy" IMF (blue) favored by 
iFardal et al. I i2007l) for reconciling the cosmic star formation history, the 
present-day ii'-band luminosity density, and extragalactic background light 
constraints. The cumulative z = evolving IMF is similar to the paunchy 
IMF, showing that it would go well towards reconciling fossil light data 
with the integrated cosmic SFH (see 0. 
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Figure 5 . Cos mic star formation history. Data points are from 
iHopkins I t2004); b est-fit relation to an updated compilation from 
Hopkins & Beac omi j20O6l) is shown as the green line. A schematic rep- 
resentation of the cosmic star formation history assuming an evolving 
IMF with 7 = 2 is sh own as the blue line, which is simply the 
Hopkins & Beacomi j2006l) fit divided by the ratio of high-mass star for- 
mation to stellar mass growth for an evolving IMF versus a constant one, 
namely the solid (7 = 2) line in Figure |3] Magenta and red points show 
the result of differentiatin g the cosmic stellar mass den sity measurements 
as a f unction redshift by iPerez-Gonzalez et al. and lElsner et al. I 

respectively. Cyan dashed line shows the fit to cosmic star for- 
m ation rate inferred from a compilation of cosmic stellar mass densities 
bv lWiUdns, Trentham & Hopkins I i2007l) . For ease of comparison, all data 
here are shown here assuming a Salpeter IMF. The evolving IMF case is in 
good agreement with direct measures of stellar mass growth. 



an d that (unlike e.g. t he merger-induced top-heavy IMF suggested 
bv lLacev et al. l2007h it appUes to the bulk of star forming galaxies 
at any epoch. 



7 IMPLICATIONS OF AN EVOLVING IMF 

What are the observational implications of such an evolving IMF 
at early epochs? Surprisingly few, as it turns out, so long as 
the massive end of the IMF remains unaffected. This is be- 
cause, as mentioned before, virtually all measures of high-2 star 
formation, be they UV, IR or radio, trace light predominantly 
from high-mass stars. Feedback energetics and metallicities like- 
wise reflect predominantly massive star output. Hence so long 
as an IMF preserves the same high mass star formation rate, 
it is expected to broadly preserve the successes of understand- 
ing the r elationship between the rest-UV properties of galax- 



les (e.g. |d ave. Finlator & Oppenheimer] I2OO6I: iBouwens et al. I 
2007h. cosmic metal pollutio n (e.g. lOppenheimer &Davell200a 



Dave & Oppenheimer I l2007h . feedback strength, and Type 11 su- 
pernova rates. 

The main difference caused by such IMF evolution is in stellar 
mass accumulation rates. High-z measures of s tellar mass ar e there- 
fore critical for testing this type of scenario. iDrorv et al. I feOOSf) 
determined the stellar mass function evolution from a K-band sur- 
vey out to 2 ~ 5, but at even moderately high redshifts the rest- 
frame light is actually quite blue, so their stellar mass estimates 
implicitly involves a signi ficant IMF assump tion. The high- 2; stellar 
mass-metallicity relation terb et al. l2006ij) has the potential to test 
the IMF; if the evolving IMF is correct, the good agreement found 



versus recent models ('Finlator & Dave I I2007I ; Isrooks et alTlbOOTl ; 
Kobavashi, Springel & White 2007) would perhaps not stand, and 
in particular many more metals would have to be driven out of 
galaxies of a given M « via outflows (i.e. the "missin g metals" prob- 
lem would get worse; IPave & Oppenheimer l2007h . Type la super- 
nova rates at high-z would be lowered; the apparent turn-down in 
Type la rates at 2; > 1 argued by Hopkins & Beacom 1 1 |2006!) as be- 
ing indicative of a long Type la delay time of 3 Gyr, may instead be 
accomodated via a more canonical short delay time coupled with 
an evolving IMF. These are examples of the types of observations 
that could potentially constrain IMF evolution. 

The star formation rates inferred for high-redshift galaxies by 
assuming present-day IMFs would have to be revised downwards 
for an evolving IMF. A quantitative estimate of the reduction in 
star formation rates from the standard to evolving IMF case is pro- 
vided by the ratio /ov//st as a function of redshift from Figure [3] 
Specifically, the "true" cosmic star formation rate (assuming the 
evolving IMF is correct) is the one currently inferred assuming a 
non-evolving IMF, divided by the curve for a particular value of 7, 
say 7 = 2. 

Figure [5] shows the cosmic sta r form ation history with 
data points compiled by iHopkins] j2004h . and the fit by 
iHopki ns & Beacom I i200q) . It also sc hematically illustrates the ef- 
fect of an evolving IMF by dividing the Hopkins & Beacom ( 200^ 
fit by the fcv/fet curve for 7 = 2. Note that Hopkins & Beacorn] 
(2006) used a Salpeter IMF to infer star formation rates, so these 
correction factors (derived for a Kroupa IMF) are not exact, but 
should be reasonably close. Figure |5] shows that the cosmic SFR 
peak shifts to later epochs in the evolving IMF case. The total stel- 
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lar mass formed in the Universe is reduced by 55% by z = and 
75% by 2 = 2, compared to a standard IMF. 

Interestingly, some current observations seem to prefer an 
IMF at early epochs that is more top-heavy or bottom-light, and 
perhap s a later epoch of cosmic sta r formation than generally be- 
lieved. iHopkins & Beacorn ] l l2006h noted that the integral of the 
cosmic star formation history exceeds the stellar m ass currently 
detec ted by a factor of two, assuming a Salpeter IMF. lBorch et al. I 
ilOOS) also noted such tension, but suggested that switching from 
Salpeter to a Chabrier or Kroupa IMF might go significantly to- 
wards reconciling the difference; little evidence exists for IMF evo- 
lution to z < 1 (Bell et al. 2007). But to higher redshifts, obser- 
vations bv lRudnick et al. I ( 12006) of stellar mass density evolution 
to 2 ~ 3 suggest that mass buildup is peaked towards significantly 
later epochs than in current mo dels that broadly match the obs erved 
cosmic star formation history. Baldrv & Glazebrook (20p 3[). up - 
dating an analysis done by Madau, Pozzetti, & Dickinson | (199g), 
found that the cosmic SFH and present-day cosmic luminosity den- 
sity are best reconciled using an IMF slope that is slightly more top- 
heavy than Salp eter (—1.15 vs. —1.35), and signifi c antly more top- 
heavy than e.g. iKroupa. Tout, & Gilmore I ( Il993h . [Heavens et al. " 
( |2004|) suggests from an archaeological analysis of SDSS spectra 
that the cosmic star format ion rate peaked at z ^ 0.6 rather than 
the more canonical z ~ 2. 1 Van Dokkum & van derMarelllioOTl) 
find that it is easier to understand early-t}^e fundamental plane evo- 
lut ion at z < 1 with a slightly top-heavy IMF (as also argued 
bv lRenzim|[2005h . Intracluster medium metallicities indicate that 
stars formed in clusters at early epochs either ha d higher yields 
than t oday, or else formed w ith a top-heavy IMF ^Portinari et al. I 
|2004) . iLucatello et aTl j2005l) and lTumlinson] i20m suggest that 
the abundance of carbon-enhanced metal-poor stars is indicative of 
a more top-heavy IMF at early times in our Galaxy. iCharv U2007h 
finds that an IMF slope more top-heavy than Salpeter is required 
to produce enough early photons to reionize the Universe. While it 
should be pointed out that a large number of observations are con- 
sistent with a universal IMF, whenever a discrepancy exists, it is 
invariably in the sense of favoring an IMF with more massive stars 
at early times. 

iFardal et al. I ilOOH) did perhaps the most careful work 
on trying to reconcile fossil light measures with redshift- 
integrated measures of stell ar m ass. Extending the analyses of 
lBaldrv& Glazebrook I i2003h and 'Hopkins & Beacom I ( |2006|) by 
including extragalactic background light (EBL) constraints, they 
found that no standard IMF (Kroupa, Chabrier, Salpeter, etc.), is 
able to reconcile at the ~ 3cr level the cosmic star formation his- 
tory, the present-day A"-band luminosity density, and the EBL. 
They examined how the IMF might be altered to make these ob- 
servations self-consistent, and found that what was required was an 
IMF that had an excess at intermediate masses (~ few i\/0) over 
standard IMFs. This IMF need not be universal, but instead could 
be a cumulative IMF of all stars that have formed in the Universe, 
since it was obtained from fossil light considerations. An example 
of a concordant IMF is their "paunchy" IMF, shown as the blue 
dot-dashed line in Figure|4l 

Using the evolving Kroupa IMF, it is possible to compute a 
cumulative IMF of all star formation in the Universe. This is done 
by integrating the evolving IMF over cosmic time weighted by the 
cosmic star formation history, which is taken to be the evolving 
IMF one (blue line) from Figure |5] The result is shown as the 
do tted cyan lin e in Fi gure |4l It is similar to the paunchy IMF 
of iFardal et al. I ( |2007|) . showing elevated formation rates of in- 
termediate mass stars, though it is slightly more top-heavy. This 



comparison is preliminary because the EBL and current if-band 
light trace stellar output differently than the cosmic star forma- 
tion history. Still, it is likely that this evolving IMF will go sig- 
ni ficantly towards rec onciling the various observations considered 
inlFardal et"an (|2007[). 



Van Dokkum ] j2007l) recently argued for a more top-heavy 



IMF at early times from the evolution of colors and mass-to-light 
ratios of early-type galaxies in clusters from z ~ 0.8 ^ 0. Pa- 
rameterizing the evolution in terms of the characteristic mass of a 
Chabrier IMF, he favored a characteristic mass that evolves by x 20 
from z ~ 4 until today; this is consistent with Mimf oc (1 + z)^ as 
proposed here. Note that his analysis only applies to stars formed 
in early-type cluster galaxies, so is not necessarily applicable to 
all galaxies. Making the assumption that it is widely applicable, 
I Van Dokkum I bOOTD determined a cosmic stellar mass growth rate 
that is broadly consistent with the evolving IMF case shown in Fig- 
ure [5] (though in detail it is more peaked towards lower z; see his 
Fig. 13). While many uncertainties are present in the his analysis, 
it is reassuring that similar IMF evolution is inferred from a com- 
pletely independent line of argument. 

The ultimate test of an evolving IMF is to directly 
measure the stellar mass buil dup in the Univers e. Recently, 
IPerez-Gonzalez et al. I ( |2007|) and lElsner et al. I ilOOH) did just that 
using deep Spitzer/IRAC observations that probe near-IR light all 
the way out to z ~ 3+. Their data of cosmic stellar mass mass 
growth rates as a function of redshift are shown as the magenta 
and red data points in Figure |5] respectively (assuming a Salpeter 
IMF). These data points were obtained by differentiating the stel- 
lar mass densities as a function of redshift published in those pa- 
pers. Although those two samples are independent, both show that 
the stellar mass growth rate is significantly different than that in- 
ferred fr om the observed cosmic star fo rmation history, particularly 
at z > 2. |Perez-Gonzalez et al. I j2007t) notes that this may favor an 
IMF weighted more towards massive stars at high redshifts. They 
mention that utilizing a Chabrier or Kroupa IMF will lessen the 
discrepancy, but such a reduction is essentially by a fixed factor at 
all redshifts (modulo minor stellar evolution differences), so even 
if the normalization is adjusted to agree with the cosmic SFH at 
z = 0, it would still yield significantly lower SFRs at z > 2. 

A major uncertainty in this comparison is that 
iPerez-Gonzalez et al. ] ( |2007|). in order to integrate up the to- 
tal stellar mass, assumed a lschechterl ( Il976 ) faint-end slope that 
is fixed at Q = —1.2, even though they cannot directly constrain 
Q at z > 2. They justify this assumption by noting that there 
is no evidence for a change in a from z ~ ^ 2. But if a 
becomes progressive ly steeper with redshift, as suggested by 
iFontana et al. I ( l2006h . it could also reconcile these observations. 
Note that it would have to steepen quickly, e.g. to a « —1.8 at 
z --^ 3 to reconcile the ~ x 3 difference, significan tly faster than 
da/dz = -0.082 o bserved bv lFontana etal. l ( |2006|) out to z ~ 4. 
lElsner et all ( l2007h also found that a does increase with redshifts 
slightly in one form of their analysis, though again not by enough 
to explain the discrepancy. Still, the difficulty of constraining the 
faint-end slope at high-z, along with inherent uncertainties in 
SED fitting, mean that these observations cannot be considered 
definitive proof of a varying IMF. But at least at face value the 
agreement with the evolving IMF case is quite good. 

Wilkins, Trentham & Hopkins (2007) independently deter- 
mined the cosmic stellar mass growth rate by compiling observa- 
tions of stellar mass densities from the literature for z ~ — 4 (note 
that they do not include the Perez-Gonzalez et al. or Eisner et al. 
data). Carefully accounting for various different systematics, they 
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fit a ICole et al. I ( 1200 ih form to the implied cosmic star forma- 
tion history (including effects of stellar mass loss), and determine 
p* = (0.014 + 0.11z)/l[l + (^/1.4)]^■^ in Moyr'^Mpc^^ This 
relation is shown as the cyan dashed line in Figure [5l it has been 
adjusted by +0.22 dex to change from the IMF they assumed (sim- 
ilar to Chabrier) to a Salpeter IMF. It agrees well with other obser- 
vations, and lies significantly below the star formation rate inferred 
from high-mass SFR indicators at high-z:. Wilkins et al. suggest that 
an evolving IMF may be the cause, and propose their own form of 
IMF evolution to explain it. As seen in Figure|5] the agreement with 
the evolving IMF proposed here is quite good. 

In short, an evolving Rroupa IMF does not obviously intro- 
duce any large inconsistencies with high-z galaxy data, and more- 
over may be preferred based on low-redshift fossil light consider- 
ations and observed cosmic stellar mass growth rates. Hence the 
amplitude evolution of the A/, —SFR relation may be yet another 
example of an observation of stellar mass buildup that suggests that 
the IMF is weighted towards more massive stars at higher redshifts. 



8 SUMMARY 

Implications of the observed stellar mass-star formation rate corre- 
lation are investigated in the context of current theories for stellar 
mass assem bly. The key point, found here and poin ted out in previ- 
ous studies jPaddi et allbOOTi : lElbaz et al.]|2007h . is that the am- 
plitude of the M, — SFR relation evolves much more rapidly since 
2 ~ 2 in observations that in current galaxy formation models. 
It is shown here that this is true of both hydrodynamic simula- 
tions and semi-analytic models, and is broadly independent of feed- 
back parameters. In contrast, the slope and scatter of the observed 
A/* —SFR relation are in good agreement with models. 

The tight A/*— SFR relation with a slope near unity pre- 
dicted in models is a generic result owing to the dominance of 
cold mode accretion, particularly in early galaxies, which pro- 
duces rapid, smooth and relatively steady infall. The slow ampli- 
tude evolution arises because star formation starts at z > 6 for 
moderately massive star-forming galaxies and continues at a fairly 
constant or mildly declining level to low-z. The large discrep- 
ancy in amplitude evolution when compared to observations from 
2 ~ 2 ^ hints at some underlying problem either with the mod- 
els or with the interpretation of data. A convenient parameteriza- 
tion of the problem is through the star formation activity parameter, 
Oaf = (A/*/SFR)/(fH - 1 Gyr), where tn is the Hubble time. 
A low value corresponds to a starbursting system, a high value to 
a passive system, and a value near unity to system forming stars 
constantly for nearly a Hubble time. In models, Osf remains con- 
stant around unity from 2 = — 4, whereas in observations it rises 
steadily from < 0.2 at 2 '--^ 2 to close to unity at 2 ~ 0. 

Several ad hoc modifications to the theoretical picture of stel- 
lar mass assembly are considered in order to match agf evolution, 
but each one is found to be in conflict with other observations of 
high-redshift galaxies. Bursts seem unlikely given the low scatter 
in Aft— SFR. Delaying galaxy formation to match a^i results in 
such a low redshift for the onset of star formation, it is in conflict 
with observations of star forming galaxies at earlier epochs. Hid- 
ing a large population of galaxies as passive runs into difficulty 
when compared to direct observations of the passive galaxy frac- 
tion in mass-selected samples at 2 ~ 2. Having an exponentially 
growing phase of star formation or "staged" galaxy formation are 
perhaps the most plausible solutions from an observational view- 
point, but it is difficult to understand theoretically how such sce- 



narios can arise within hierarchical structure formation. Hence if 
observed star-forming galaxies are quiescently forming stars and 
represent the majority of galaxies at 2 ~ 2, as other observations 
seem to suggest, then it is not easy to accomodate the low Qsf in- 
ferred from the A/, —SFR amplitude. 

Systematic uncertainties in Al, or SFR determinations could 
bias results progressively more at high-2 in order to mimic evolu- 
tion in Qsf. Various currently debated sources are considered, such 
as the contribution to near-IR light from TP-AGB stars, extinc- 
tion corrections, assumptions about star formation history, AGN 
contamination, and PAH emission calibration. It may be possible 
to concoct scenarios whereby several of these effects combine to 
mimic a^f evolution, but there are no suggestions from local obser- 
vations that such scenarios are to be expected. Hence if systematic 
effects are to explain the low Qsf at high-2, it would imply signifi- 
cant and unexpected changes in tracers of star formation and stellar 
mass between now and high redshifts. 

A solution is proposed that the stellar initial mass function be- 
comes increasingly bottom-light to higher redshifts. Several lines 
of arguments are presented that vaguely or circumstantially favor 
such evolution, though no smoking gun signatures are currently 
known. A simpl e model of IM F evolution is constructed, based 
on the ansatz bv lLarsonI j 19981) that the minimum temperature of 
molecular clouds is reflected in the characteristic mass of star for- 
mation (AfiMp) where the mass contribution per logarithmic mass 
bin is maximized. The minimum temperature may increase with 
redshift owing to a hotter cosmic microwave background, more 
vigorous star formation activity, or lower metallicities within early 
galactic ISMs. In order to reconcile the observed asi evolution with 
theoretical expectations of no evolution, an evolving IMF of the 
form 

'^^ (X M-°'^ for M < MiMF, 



dlog Af 
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for M > AfiMp; 



= 0.5{l + zfMQ 



is proposed. The exponent of AfiMP evolution is constrained by 
requiring no Qsf evolution, through careful modeling with the PE- 
GASE.2 population synthesis code. While the exact form of IMF 
evolution is not well constrained by present observations, what is 
required is that the IMF has progressively more high-mass stars 
compared to low-mass at earlier epochs, and that this IMF applies 
to the majority of star forming galaxies at any epoch. It is worth not- 
ing that this evolving IMF is only constrained out to 2 ~ 2 from the 
Af, — SFR relation, though it yields predictions that are consistent 
with other observations out to 2 ~ 4. Extrapolating such evolution 
to higher redshifts is dangerous, since no observational constraints 
exist and the precise cause of IMF evolution is not understood. 

Implications of such an evolving IMF are investigated. By 
leaving the high-mass end of the IMF unchanged, recent successes 
in understanding the connections between high-mass star forma- 
tion, feedback, and metal enrichment are broadly preserved. The 
cosmic stellar mass accumulation rate would be altered compared 
to what is inferred from cosmic star formation history measure- 
ments using a standard IMF. It is shown that an evolving IMF is at 
face value in better agreement with direct measures of cosmic stel- 
lar mass assembly ( |P erez-Gonz alez et al. 2007; Eisner et al. 200^; 
Iwilkins. Trentham & Hopkins I boOTi) . Furthermore, the evolving 
IMF goes towards relieving the generic tension between present- 
day fossil light measures versus observations of the cosmic star 
formation history. In particular, the paunchy IMF favored by 
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iFardal et al.] ( l2007h in order to reconcile the observed cosmic star 
formation history, present-day K-hznd luminosity density, and ex- 
tragalactic background light constraints, is qualitatively similar to 
the cumulative IMF of all stars formed by today in the evolving 
IMF case. Individually, each argument has sufficient uncertainties 
to cast doubt on whether a radical solution such as an evolving 
IMF is necessary. But taken together, the il/, — SFR relation adds 
to a growing body of circumstantial evidence that the ratio of high- 
mass to low-mass stars formed is higher at earlier epochs. 

It is by no means clear that IMF variations are the only vi- 
able solution to the M, — SFR dilemma. The claim here is only that 
an evolving IMF is an equally (un)likely solution as invoking un- 
known systematic effects or carefully crafted star formation histo- 
ries in order to explain A/, —SFR evolution. It is hoped that this 
work will spur further efforts, both observational and theoretical, 
to investigate this important issue. 

Future plans include performing a more careful analysis of 
the evolving IMF in terms of observable quantities, in order to 
accurately quantify the impact of such IMF evolution on the in- 
terpretation of UV, near-IR, and mid-IR light. Also, it is feasible 
to incorporate such an evolving IMF directly into simulation runs 
to properly account for gas recycling in such a scenario. Finally, 
including some feedback mechanism to truncate star formation in 
massive systems may impact the A/,— SFR relation in some way, 
so this will be incorporated into the hydro simulations. 

Observationally, pushing SFR and A/* determinations to 
higher redshifts is key; a continued drop in Osi to z ^ 3 would 
rapidly solidify the discrepancies with current models, and would 
also generate stronger conflicts with other observations of high-z 
galaxies. Assessing the AGN contribution and extinction uncertain- 
ties from high-z systems is critical for accurately quantifying the 
light from high-mass star formation, for instance through the use 
of more direct star formation indicators such as Paschen-a. Push- 
ing observations further into the mid-IR such as to 70/i, past the 
PAH bands at z ~ 2, would mitigate PAH calibration uncertainties 
in SFR estimates. Obtaining a large sam ple of spectra fo r typical 
high-2 star-forming systems (like cB58: iPettini et al.ni2000) would 
more accurately constrain the SED than broad-band data. All of 
these programs push current technological capabilities to their lim- 
its and perhaps beyond, but are being planned as facilities continue 
their rapid improvement. 

In summary, owing to the robust form of star formation histo- 
ries in current galaxy formation models, the Af«-SFR relation rep- 
resents a key test of our understanding of stellar mass assembly. 
Current models reproduce the observed slope and scatter remark- 
ably well, but broadly fail this test in terms of amplitude evolution. 
Whether this reflects some fundamental lack of physical insight, or 
else some missing ingredient such as an evolving IMF, is an issue 
whose resolution will have a significant impact on our understand- 
ing of galaxy formation. 
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